Ionizing radiation from hydrogen recombination strongly suppresses the hthium 

scattering signature in the CMB 
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It has been suggested that secondary CMB anisotropies generated by neutral lithium could open a 
new observational window into the universe around the redshift z ~ 400, and permit a determination 
of the primordial lithium abundance. The effect is due to resonant scattering in the allowed Li I 
doublet (2s - 2p 

^Pi/2,3/2), so its observability depends on the formation history of neutral 
lithium. Here we show that the ultraviolet photons produced during hydrogen recombination are 
sufficient to keep lithium in the Li 11 ionization stage in the relevant redshift range and suppress 
the neutral fraction by ~ 3 orders of magnitude from previous calculations, making the lithium 
signature unobservable. 

PACS numbers: 98.62.Ra, 98.70.Vc, 26.35.+C, 98.80. Es 



I. INTRODUCTION 

There are few observable traces of structure in the 
universe between the cosmic recombination at redshift 
z « 1100 and the redshifts z < 7 accessible to quasar ab- 
sorption line measurements Currently our main ob- 
servational constraints on this redshift range come from 
the early integrated Sachs- Wolfe (ISW) effect Pj and the 
large-scale Thomson scattering features in the CMB 3] . 
These may be supplemented in the future by measure- 
ments of the Hi 21 cm hyperfine line 0, S IE S 
However, none of these probes can provide detailed in- 
formation in the redshift range 200 < z < 500: the free 
electron abundance is too small to produce significant 
Thomson scattering, the gas is in thermal equilibrium 
with the CMB so there is no 21 cm signal, and the early 
ISW effect has a blackbody frequency dependence that 
has no redshift information, making the detection of the 
effect model-dependent. 

The series of papers 0, 0, 0| introduces the prospect 
of using anisotropies in the CMB generated by scattering 
in the resonance lines of neutral lithium to study large 
scale structure in the redshift range z w 200-500. Be- 
cause it involves a resonance, the frequency dependence 
of the Li I scattering feature can be mapped into a red- 
shift dependence just as is done for the Lyman-a forest 
and is proposed for high-redshift 21 cm studies. The fre- 
quency dependence deviates strongly from a blackbody, 
and thus should be distinguishable from the CMB. This 
idea is also attractive because it provides a potentially 
clean method of determining the abundance of primor- 
dial lithium, before it has been destroyed by nuclear pro- 
cesses in stars. A precision measurement of the lithium 
abundance from Big Bang Nucleosynthesis (BBN) would 
be an important consistency check of the cosmological 
model, and would constrain extensions to standard cos- 
mology such as inhomogeneous BBN 0, 0| , which 
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tend to increase the ''Li abundance. 

The primordial abundance of lithium is extremely 
small, Xu « 10-i°-10-9 but the high integrated 
cross section for excitations of the Li i doublet 2s — 
2p ^Pi/2,3/2 suggests the possibility of producing a signa- 
ture within the reach of studies by the Planck, BLAST, 
and EDGE experiments - generally experiments with 
sensitivity around the 300 fim band 9|. The Lii dou- 
blet has an energy of 1.847 eV (A = 6708 A), so that 
the resonant feature would appear to us today at A = 
0.6708(1 -I- z) ^m for scattering at redshift z. By compar- 
ison, the lowest single-photon excitation for Li 11 occurs 
at 61.3 eV|33, where there are essentially no CMB pho- 
tons to scatter (and to which the pre-reionization uni- 
verse is opaque due to hydrogen and helium photoion- 
ization). Hence, an observable lithium signature on the 
CMB relies on the formation of neutral lithium after hy- 
drogen recombination. 

Previous work on the lithium signature |E llOl ITlj em- 
phasizes that with a thermal CMB spectrum, lithium 
recombines around z = 400 to form a significant neu- 
tral fraction. The neutral lithium gas has an optical 
depth that is high enough to generate a secondary CMB 
anisotropy, visible today around 268 fim. Important to 
the result is that there are too few ionizing photons in 
a purely thermal CMB to keep lithium ionized. A sepa- 
rate line of work discusses the non-thermal radiation field 
generated by hydrogen recombination IT^ Il8j , which 
contains ultraviolet (UV) photons from both Lyman-a 
production and 2s — > Is two-photon decays. These pho- 
tons have a maximum energy of 10.1988 eV, nearly dou- 
ble the Li I ionization threshold, XLi = 5.3917191 eV 
Here, we show that this non-thermal ionizing radiation 
is sufficient to keep lithium almost completely ionized up 
to late times, and to hold the optical depth through the 
resonance below 10^** in the standard cosmology. 

In this paper, we first calculate the spectral distortion 
in Sec.ni using the three-level approximation for hydro- 
gen recombination. Sec. Illll then discusses the photoion- 
ization and recombination rates of lithium in the presence 
of the spectral distortion. In Sec. IIVI we solve the rate 
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equation for lithium's ionization history, and argue that 
it remains mostly ionized. Sec.lVloutlines the significance 
this has for a lithium CMB signature. Finally, because 
this is a non-standard use of the three-level approxima- 
tion, we justify it subject to some possible concerns in 
Sec. IVII and sum up the main results in section rVIII 



II. ATOMIC HYDROGEN SPECTRAL 
DISTORTIONS 

Neutral lithium can not form if the radiation field con- 
tains a sufficient number of photons above the ioniza- 
tion threshold. This makes photoionization dependent 
on the details of high energy tail of the photon distri- 
bution {E > XLi), and the photoionization cross section 
above threshold. The cosmological radiation has nearly 
a blackbody distribution '20'|, except for a small popu- 
lation of non-thermal distortion photons from hydrogen 
recombination in the high energy tail of the distribution 

[Ulillil. 

The most important sources of non-thermal radiation 
during hydrogen recombination are the Hi two-photon 
transition (2s — > Is) and the Hi Lyman-a transition 
(2p Is). Very few distortion photons are produced 
in the higher-order Lyman lines because such photons 
are efficiently destroyed by absorption followed by decay 
to an excited state of Hi. Also, the Balmer, Paschen, 
etc. lines of H i, as well as nearly all of the photons pro- 
duced by direct recombination to H i excited states have 
energies below 5.39 eV, and so cannot ionize Lii. The 
populations of the 2s and 2p states are expected to vary 
out of equilibrium at late times. The bearing this has on 
the formation history of neutral lithium is discussed in 
Sec. EH 

We will see that the non-thermal photons from n = 
2 — > 1 decays of H i have an abundance of 10^^ relative to 
the the peak of the thermal distribution, but exceed the 
CMB photon population above the ionization threshold 
of lithium 117)], on the high-energy tail of the CMB. The 
ionizing radiation from these decays must be included in 
a calculation of lithium's ionization history. 

A natural quantity to track is the number of photons 
per hydrogen, per log-energy. Written in terms of the 
phase space density J\f and specific intensity Jjy, 
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where nn is the hydrogen number density and h and c 
are Planck's constant and the speed of light, respectively. 
This choice also has the advantage that it is conserved 
along a trajectory, unlike J^. The photon number densi- 
ties and baryon number densities both dilute by (1 -I- z)'^ . 
With these units, the high-energy Wien tail of the black- 
body distribution scales as, 
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(After z « 200, matter and radiation are decoupled, so 
their thermal histories depart and must be tracked sep- 
arately as a radiation temperature Tr, and matter tem- 
perature Tm-) From the integral over the asymptotic 
distribution, there are ~ 80 photons capable of ioniz- 
ing lithium, per lithium atom at z — 500, and only 
« 1 X 10~^^ at 2; = 300 from the thermal distribution. 
Purely thermal CMB radiation is not enough to keep 
lithium ionized after z = 400, but the distortion extends 
to much higher energies. 

We find the distortion by calculating the rates of two- 
photon and Lyman-a decays, and transport the spec- 
trum of photons that is generated from reactions of these 
rates and the known emission profiles. The method works 
post-facto because it takes the ionization history of hy- 
drogen and its derivatives, and then calculates the radi- 
ation field. Ideally, the evolution of the radiation field 
would be included in the recombination calculation, self- 
consistently, rather than computed later. Results from 
a recombination calculation with a self-consistent radia- 
tion field are expected to be different at the percent level 
because of feedback and other effects . Here, we want 
to determine whether lithium remains ionized, and do 
not need anything so nuanced. 

Assuming that only the two n = 2 decays of hydro- 
gen to its ground state contribute to the distortion, all 
the rate information can be extracted from the total hy- 
drogen recombination rate in the three-level approxima- 
tion filEil, 
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where Xe is the fraction of electrons relative to hydrogen, 
an is the case B recombination rate of hydrogen, (3h is 
the effective case B photoionization rate, 
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Aa-y is the 2s 



Is two-photon decay rate, X2s the ioniza- 
tion energy from the 2s state, and K = A|j^^/[87r7?(z)] 
is the Peebles inverse resonance escape rate from the 
Lyman-a line. The form of the three-level rate equation 
is suggestive: the term on the second line is the fraction of 
excited hydrogen atoms that make the transition 2p — > Is 
or 2s —> Is, relative to the atoms that are photoionized. 
The term in the numerator displays a two-photon piece 
and a Lyman-a piece, while the denominator also has a 
contribution from photoionization. This makes it easy to 
decouple the relative contributions of each process to the 
overall hydrogen recombination rate: call the fractional 
rate of change in Xe during hydrogen recombination due 
to two-photon and Lyman-a transitions {dxe/ dz)\2-y and 
(dXe/dz)\Lya, respectively. 

The solutions to the three-level rate equation and 
evolution of the matter temperature are evaluated by 
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reef ast, developed in Seager et al. |2l], |23- We have 
modified it to return the fractional rates, {dxe/ dz)\2j 
and {dxe/ dz)\Lya as a function of z. Throughout, the 
calculations are based on the cosmology Q.m — 0.27, 
= 0.046, Hq = 70 km s"! Mpc~\ Tcmb = 2.728 K, 
and Yhc = 0.24. 

To understand how photons generated during hydro- 
gen recombination influence the ionization history of 
lithium, they must be transported to later redshifts. 
Since Thomson scattering does not significantly change 
the photon energy at the redshifts of interest [Sg, trans- 
port amounts to redshifting the emitted distributions. 
The set of distortion photons around some energy E is 
the sum over all photons which were produced at an ear- 
lier time (zcm) with energy i?om which could have red- 
shifted to the energy E by redshift z. This defines an 
integral along the line of sight of the contribution of a 
particular process to the radiation field. 



Evolution of the radiation field 



' process 



(5) 

where 0(i?om) is the emission profile of the process, eval- 
uated at the energy of the radiation at redshift Zom that 
will have redshifted down to an energy E by redshift z, 
E{z)cai = £'(1 -I- Zom)/(l + z) . The combination of the 
profile and frequency suggests a new quantity: 
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The integral Eq. Q along a line of sight is equivalent to 
an integral over energy under the change of variables: 

i?e,n = E^l±^ ^ dz,^, = (1 + Zem)^. (7) 
(1 + Z) Ecm 

The integral along the line of sight adds contributions 
from increasingly higher redshifts, with higher emission 
energies Ecm- Because the emission profile is the same 
as time goes on, this is equivalent to adding up different 
contributions over the profile evaluated at progressively 
higher energies. The change of variables amounts to con- 
verting this indirect way of summing over redshifts into 
a direct, convenient sum of the profile and rates over en- 
ergy. 
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Here, Zem is the redshift at which light was emitted with 
energy E^m that will redshift to energy E at some later 
redshift z, 
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We use the polynomial fit to the Spitzer & Greenstein |23 | 
two-photon profile provided in 24], and assume a delta 
function for the Lyman-a profile. Results from the cal- 
culation are shown in Fig. ^ 
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FIG. 1: The distortion field generated by n = 2 decays of 
hydrogen during recombination dwarfs the high-energy tail of 
the CMB (shown in black dotted line). The ionization energy 
of hthium is shown to emphasize that even at late times, there 
is a large number of distortion photons above the ionization 
energy. 



III. LITHIUM RECOMBINATION AND 
PHOTOIONIZATION 



The photoionization rate of Li i has a thermal contri- 
bution from the CMB and an non-thermal contribution 
from the integral of the cross section over the distortion 
field, /^Li = /3cMB + /3dist ■ The thermal photoionization 
rate is related by the Milne relations to the recombina- 
tion rate aLi, 



/3cMB 



5LiI 



aLi(r,.)exp 



XLiI 



(27rmeA:T,.)3/2 



(10) 

The Milne relations assume Boltzmann-distributed ex- 
cited states at temperature T^. Yet, the distortion pho- 
tons may increase the population of excited states, in- 
creasing the ionization rate from CMB photons. If this 
effect is important, it can only increase /3cmBj and thus 
reduce the neutral lithium abundance. Ionization by a 
distortion photon occurs predominantly from the ground 
state. For the case A Liii— >Lii recombination rate, aLi, 
we use a fit from Verner and Ferland (25(1, 




where the coefficients are, a — 1.036 x 10 i^'-^'^o ^ 



To\ y To) y y Ti 

(11) 

cm"' s 

b = 0.3880, To = 1.077 x 10^ K, and Ti = 1.177 x 10^ K. 
The fit is accurate to better than 4% (with maximum 
error less than 6%) [l^ between 3 K and 10^ K. 

The non-thermal contribution to lithium photoioniza- 
tion from photons emitted during hydrogen recombina- 
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FIG. 2: The photoionization rate of lithium as a function 
of redshift, showing the contribution of non-thermal (distor- 
tion) photons, beginning to dominate shortly after hydrogen 
recombination. The non-thermal rate is fast compared to the 
inverse Hubble time, allowing a steady-state solution. 

tion is the integral over the cross section, 

r^^y dE' f \ 
/3dist = y —\r2^{E')^ri^y^{E')Y^,(E')c. (12) 

We use a fit to the photoionization cross section for 
Lil from that has been smoothed over resonances. 
(There are no resonances in our energy range, so this is 
not a concern.) The fit is based on low-energy cross sec- 
tions from the Opacity Project [23| and high-energy cross 
sections from a Hartree-Dirac-Slater method (though the 
high-energy cross section is unimportant here because 
there are no photons above Ej^ya = 10.2 eV): 

a{y = E/Eo) = fTo(2/-l)22/°-^^-^-^(l + y^)-^, (13) 

where ctq = 1.87723 x lO"*' cm^ s'^/c, P = 4.895, and 
Ha — 15.01. The results for the non-thermal and thermal 
contributions to lithium photoionization as a function of 
redshift are shown in Fig. |21 

IV. IONIZATION HISTORY AND SCALING 

Lithium chemistry and the interaction between lithium 
ions, lithium hydride, etc. can be very complicated .28j. 
but we will follow Stancil et al. Q in tracking only pho- 
toionization and radiative recombination, the two most 
dominant reactions, 

Li+ +c" ^ Li + 7. (14) 

The next leading order reaction is radiative charge trans- 
ferif2l, 

Li + H+ ^ Li+ +H + 7, (15) 



FIG. 3: Lithium recombination history in the presence of an 
ionizing flux from the n = 2 hydrogen recombinations to the 
ground state. Lithium recombines much more slowly than for 
a purely thermal CMB, as (1 + z)~^^'^ , and remains largely 
ionized up until late times. This reduces the optical depth, 
and significantly reduces the imprint left on the CMB. 

which is two orders of magnitude slower than Eq. H14() 
if distortion photons are neglected. Photoionization 
from distortions will only speed up Eq. (|14|l relative to 
Eq. H15|) . making charge transfer less important. 

The rate equation describing the evolution of the neu- 
tral lithium fraction fu based on the recombination au 
rate and photoionization rate Pu is 

^ - -(l + z)-lff(z)-lKine(l-/Li)-/3Li/Li]. (16) 

The solution of the rate equation Eq. H16() is shown in 
Fig. 13^ 

Lithium recombination according to Eq. (|16|) goes 
through four major stages: 

1. A Saha stage, where the lithium ionization states 
are in thermodynamic equilibrium; 

2. A primary steady-state stage, in which the spectral 
distortion from the bulk of hydrogen recombination 
at z ~ 1300 keeps lithium ionized; 

3. A residual steady-state stage, in which the spec- 
tral distortion from z ~ 1300 has rcdshiftcd to be- 
low the Li I ionization threshold, but lithium is kept 
ionized by UV radiation from residual hydrogen re- 
combinations; and 

4. A freeze-out stage, in which the time scale for all 
relevant reactions becomes larger than the Hubble 
time, and /li approaches a constant. It is dur- 
ing this stage that the steady-state approximation, 
Eq. ((TH|) . breaks down. 

In the real universe, the first stars eventually turn on, 
increasing the UV background-probably ionizing essen- 
tially all of the intergalactic lithium. This occurs too late 
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(i.e. too low z) to be of interest for the lithium feature 
in the CMB. 

We now discuss each of the four stages in more detail. 



A. Saha stage (z > 700) 

Initially, Li i and Li ii co-exist in a Saha equilibrium, 
where many thermal ionizing photons are present and the 
distortion is negligible. Note that the neutral fraction 
during this phase is extremely small, /li < 10~^. The 
Saha equation in this case is 



fu = n. 



( 



3/2 



„X-L.i/(feT,) 



\ 2'KmekTr 

This equation is valid until z ~ 700. 



(17) 



B. Primary steady-state stage (400 < z < 700) 

Around z ^ 700, enough CMB photons redshift be- 
low the Li I ionization threshold so that a small number 
of distortion photons dominate the Li i photoionization 
rate. During this second stage, the neutral fraction con- 
tinues to increase, but it lags far behind the Saha pre- 
diction. Although thermodynamic equilibrium no longer 
applies in this redshift range, the lifetime (3^^ of a neu- 
tral lithium atom is still small compared to the Hubble 
time H^^ . Thus a "steady state" occurs in which the net 
Li I production rate dfu/dz is much smaller than either 
the recombination or photoionization rates in Eq. (|16|l . 
Thus we have 



fu ~ aune/iPu + aune) ~ aurie/P: 



Li- 



(18) 



[This steady-state solution for lithium is analogous to the 
"ionization equilibrium" of H I and H"*" in the Lyman-a 
forest. In both cases the recombination rate to form the 
neutral species (Lii or Hi) is counterbalanced by pho- 
toionization by UV photons from a non-thermal source, 
which consists of hydrogen recombination radiation in 
the case of lithium versus stars and quasars in the Lyman- 
a forest. Neither case represents a true thermodynamic 
equilibrium because the radiation field is non-thermal.] 



C. Residual steady-state stage {20 < z < 400) 



These act as a source of UV photons, keeping the neu- 
tral lithium fraction small. 

During the period of residual recombination, the 
steady-state solution for the neutral fraction can be 
simplified into a scaling law by making two main as- 
sumptions: matter temperature changes adiabatically 
(T oc a~^) after thermal decoupling at 2 ~ 200, and the 
relic electron fraction Xe is nearly constant after freeze- 
out. The low temperature limit of the recombination 
rate coefficient in Eq. Hll|l for adiabatic matter gives 
Q!Li oc T^^/^ oc a. The photoionization rate is then 



Pu ann^x^t oc a 



(19) 



(a oc t^/'^ for matter domination, H oc t^^. The number 
density squared scales as oc a~^.) Likewise for hydro- 
gen, an oi a, so that the neutral lithium fraction is given 
by, 
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(20) 

The /li oc c?!"^ scaling law is shown in Fig. O and does 
a remarkably good job of reproducing the full numerical 
solution in the range from thermal decoupling at z ~ 200 
to lithium freeze-out at z ~ 20. 



D. 



Freeze-out < 20) 



As the universe continues to expand, the number of 
residual recombinations drops. Consequently the UV 
background becomes fainter, and eventually the lifetime 
of a lithium atom exceeds the Hubble time, fi^l ^ . 
Here, the neutral fraction of lithium breaks away from 
the c?!"^ scaling and approaches a constant (see Fig.O. 

At some point the freeze-out is interrupted by cosmic 
reionization caused by the first stars and quasars. Cur- 
rent observations "3 0*1 indicate that hydrogen was reion- 
ized at 2 > 6; lithium was presumably ionized even ear- 
lier because the pre-reionization universe was transpar- 
ent to photons with energies of 5.4-13.6 eV (except in 
the Lyman lines), which are capable of ionizing lithium. 
Regardless, this redshift range is not of interest for ob- 
servations of the lithium resonance line scattering signal: 
scattering at z = 20 would be observed at a wavelength 
of 14 /im, where the CMB radiation is many orders of 
magnitude fainter than foreground emissions. 



When hydrogen recombination freezes out at z ^ 800, 
the residual ionization fraction Xe ^fewxlO^"* changes 
very little. The Lyman-a photons produced at z ~ 800 
are cosmologically redshifted and by 2 400 are no 
longer capable of ionizing Li i. Thus one might naively 
expect that after 2 ~ 400, lithium would fully recom- 
bine to Li i. As can be seen from Fig. |31 this is not the 
case: even after hydrogen recombination is "frozen out," 
there are still a few residual recombinations occurring. 



V. SIGNATURE IN THE CMB 

The observability of the lithium signature in the CMB 
is determined principally by the optical depth through 
the Li I doublet. As noted in previous work [oL IToL [T]| . 
the major effect of lithium on small angular scales 
is to suppress the CMB temperature and polarization 
anisotropics by a factor of exp(— tli), where tli is the 
optical depth through the Lii doublet. This suppression 
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factor depends on the scattering redshift Zg , which in turn 
depends on the wavelength of the observation through 
the usual relation Zg = Aobs/(6708A) — 1. There are also 
intermediate-scale effects ^ 100), namely the Doppler 
temperature anisotropy due to motion of the gas at the 
lithium scattering redshift, and the polarization induced 
by Li I scattering. In the optically thin regime tli <C 1 
of interest here, the magnitude of the effect as measured 
by, e.g. the change ACf in the CMB power spectrum, is 
proportional to tli. This section is primarily concerned 
with computing the optical depth tli. We will then argue 
by scaling from the results of Stancil et al. 9] that the 
signature is too small to be observed. 

The optical depth of neutral lithium in terms of the 
neutral fraction, neglecting the small fraction occupying 
the upper states is ^ 2lL iSL] 



TLi(Aobs 
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(21) 
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FIG. 4: A comparison of the optical depth through the reso- 
nance of neutral lithium, with and without the ionizing flux 
from hydrogen recombination. 



Here the sum extends over all excited levels accessible 
from the ground level, and the scattering redshifts are 
given by Zs^i — Aobs/-^ii ~ 1- Among those, only the 
2p J — 1/2 and J — 3/2 excitations at 1.8478 eV con- 
tribute for photons at low frequencies (the next allowed 
excitation is to 3p ^Pi/2,3/2 at 3.7 eV, where dust compli- 
cates observations), and each has an Einstein coefficient 
of 3.72 X 10'' s-\ wavelengths of X1/2 = 6707.76 A and 



A 



3/2 



— 6707.91 A, and degeneracies gi/2 — 2 and 



The degeneracy of the ground level 



.93/2 = 4 n Ei. 

is 51 = 2. 

We also need to know the lithium abundance in 
order to compute the optical depth. There is some dis- 
cord about the value of Xu ■ Recent standard BBN abun- 
dance calculations based on the R-matrix method for low 
energy nuclear cross sections yield Xu = i^-^^-o'.^) x 
10~^° |32, while observations of low-metallicity stars 
give Xu = (1.23l^j^) x 10"!° |3|. The BBN value 
of Xu = 4.15 X 10"-^° is used throughout this analysis. 
(If the measurement based on stars is correct, the opti- 
cal depth is correspondingly reduced, tu oc Xu-) The 
optical depth with and without the spectral distortions 
is shown in Fig. 0] 

The optical depths shown here are 3-4 orders of magni- 
tude smaller than those found by Stancil et al. Q in their 
Model A (which most closely matches the cosmology used 
here) in the interesting redshift range 200 < z < 500. In- 
stead of a 20% change in the CMB power spectrum at 
A = 268 /im {zs = 400), we have tu ^ 10~^ and find a 
fractional change of the order of \ACe/Ci\ ^ 2 x 10"^. 
A similar conclusion holds at other redshifts because we 
find that tu is always < 10^^. Thus, with the standard 
BBN abundance of lithium, the temperature or polariza- 
tion due to the lithium feature never exceeds about 10"^ 
of the primary CMB anisotropy. This is extremely diffi- 
cult to measure even in principle, and is well beyond the 
capabilities of current experiments in the relevant wave- 
length range. 



In practice, the lithium signature must also be disen- 
tangled from the far infrared background (FIB). Zaldar- 
riaga & Loeb 10] estimate that the FIB fiuctuations are 
comparable to the lithium signal in polarization and dom- 
inate over lithium in temperature anisotropy at Zg ~ 500, 
if Tu is of order unity. At lower Zg, the FIB-to-CMB ratio 
increases dramatically. Since we calculate tu < 10^**, it 
follows that the lithium signal should be at most ^ 10^^ 
of the FIB fiuctuations. Thus the lithium signal will 
be undetectable unless both the experimental capabili- 
ties improve dramatically, and the FIB fiuctuations turn 
out to be much smaller than predicted or very simple to 
model. 

It is conceivable that the lithium abundance is greater 
than that predicted by standard BBN, with some inho- 
mogeneous models giving Xu one to two orders of mag- 
nitude larger. This is only sufficient to bring tu up to 
the 10~'^-10~^ range, which is still much smaller than 
previous predictions. 



VI. SUBTLETIES 

While the three-level approximation provides a con- 
venient way to estimate the spectral distortion from hy- 
drogen recombination, it is a non-standard application of 
an approach that was developed largely to calculate the 
free electron fraction accurately. In this section, we jus- 
tify its application to the problem of lithium subject to 
concerns about lithium feedback, i.e. the possibility that 
Li I photoionizations could absorb the spectral distortion 
(Sec. I VI A|) . Hi 2s and 2p falling out of equilibrium at 
late times, influencing the production of the distortion 
and the accuracy of case B hydrogen recom- 



rsec lyiBt 

bination at late times fSec. IVI Hjl . 
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A. Lithium feedback 

As lithium recombincs to the ground state, it absorbs 
photons from the distortion field. In the most extreme 
case of this, all of the distortion photons are absorbed 
by a small number of lithiums, accelerating lithium re- 
combination. We would like to show that this "lithium 
feedback" is negligible, and that it is realistic to use the 
distortion field set up by hydrogen recombination alone. 
This can be shown by checking that the total number of 
photons absorbed by lithium that could have ionized a 
Li I atom at redshift z is much less the total number of 
distortion photons. The number of photons per H nucleus 
that are absorbed (Xrcmovod) is equal to, and bounded by. 



c 



H{z) 1 + z 

initial /Li I final } 



H{z) l + z 



1 



, (22) 



where the endpoint of the integral is the maximum red- 
shift where photon from a Lyman-a decay can be created 
so that it reaches z with enough energy to ionize lithium. 



Ze(z) 



E 



Lya 



XU 



(23) 



By comparison, the number of photons above the Li i ion- 
ization threshold in the distortion per H nucleus (Xhcm) 



-'Afield — 



j-^l '^distortion \-t^ j ^ ) 
XLi ^ 



(24) 



The two integrals, Eqs. H22|l and l|24|l . are shown in Fig.El 
The number of photons removed from the distortion is 
several orders of magnitude less than the total distortion 
population. This confirms that absorption of the distor- 
tion photons by lithium can be neglected. 



Removal of photons from tfie spectral distortion 
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FIG. 5: The maximum number of photons that lithium re- 
combinations can remove from the distortion population is 
significantly less than the total distortion population, based 
on XLi ~ 4.15 X 10"^". The distortion spectrum is well- 
approximated by just those photons which are generated by 
hydrogen, without lithium feedback. 



spectral distortion. If at late times, the majority of the 
decays are via the two-photon mechanism, then the ion- 
izing flux will not be as hard as when production is a 
mixture of two-photon and Lyman-a, or all Lyman-a. 
Here we look at the worst and best case scenarios where 
the distortion is produced by only two photon processes, 
or only Lyman-a processes, shown in Fig.|Sl This is a lit- 
mus test for whether 2s and 2p being out of equilibrium 
has any bearing on the conclusion: the lithium ioniza- 
tion history lies somewhere between both of these cases. 
As can be seen from the figure, even the two extremes 
do not alter the conclusion that lithium recombination is 
dramatically inhibited. 



C. Accuracy of recombination coefficients 



B. 2s— 2p non-equilibrium 

The second primary concern is that the spectral distor- 
tion derived from rates in the three-level approximation 
may not be accurate because it assumes that 2s and 2p 
are in equilibrium in hydrogen. The reactions that keep 
2s and 2p in equilibrium involve absorption and emis- 
sion of Balmer and Paschen line photons. These are 
expected to slow down not long after hydrogen recom- 
bination when the CMB photons redshift to sufhciently 
low energies that they cannot excite Balmer transitions. 
This means that 2s and 2p could fall out of equilib- 
rium because the ratio of effective recombination rates 
ai^J I (x^2s^ differs from the statistical ratio g2p/g2s = 3. 
This becomes a concern for how well we understand the 



The third possible concern is that the case B hydrogen 
recombination coefficients used in reef ast from Pequig- 
not et al. [S^l will not scale properly below 40 K, where 
residual hydrogen recombinations can ionize lithium (the 
coefficients are accurate at high temperatures during the 
bulk of recombination) [2^ |33 . The case A hydrogen re- 
combination rates in Verner and Ferland ^2S\ are accurate 
down to 3 K, with a maximum error of 6%. These can be 
converted to accurate case B rates that are similarly ac- 
curate down to 3 K by subtracting recombinations to the 
Is level. The fitting form for case B hydrogen recombi- 
nation of Pequignot et al. used in reef ast departs from 
this more accurate calculation by less than ~ 10% for 
temperatures above 3 K. The recombination coefficients 
used in reef ast are accurate to several percent over the 
redshift range in question. 
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The effect of far non-equilibrium distortion production 
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2s, 2p equilibrium 
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FIG. 6: For a spectral distortion produced by only Lyman-a 
decays, the spectral distortion is much harder, and lithium re- 
combines less readily at early times. Alternately, for a distor- 
tion produced by only two photon decays, the softer radiation 
allows lithium to recombine more efficiently. This simulates 
the effect of hydrogen 2p and 2s falling out of equilibrium in 
its most extreme form. In both cases, lithium remains largely 
ionized up to late times. 



VII. CONCLUSIONS 

We have shown that the non-thermal UV flux from hy- 
drogen recombination is sufheient to keep lithium almost 
completely ionized throughout the cosmic Dark Ages. 
Previous studies have only considered the thermal tail 



of the CMB, which is not enough to ionize lithium at 
times later than z « 400-500. The UV flux from hy- 
drogen recombinations dramatically reduces the optical 
depth through the Lii doublet and the prospect for ob- 
serving a lithium signal in the CMB. This is unfortunate 
because the lithium scattering signal, had it been de- 
tectable, would have been a valuable cosmological probe 
in the 200 < z < 500 redshift range and would have 
improved our understanding of BBN physics by permit- 
ting a clean measurement of the primordial lithium abun- 
dance. 

UV radiation from hydrogen recombination may also 
be important for the production of molecules such as H2, 
LiH, and HD because some of the the intermediate prod- 
ucts, catalysts and final states can be destroyed by UV 
radiation. H2 and HD have been suggested as coolants for 
gas clouds in the early universe which could be important 
for the formation of the first baryonic structures. Cool- 
ing processes will also be crucial in our understanding of 
upcoming high-redshift 21 cm surveys, both directly if 
an H I signal from minihalocs is observed 36] , and indi- 
rectly because cooling affects the formation of UV sources 
such as stars. These results highlight the importance of 
developing a code to track the radiation field during re- 
combination, and we hope to discuss a new algorithm for 
this in a future paper. 
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